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1.0  SURFACE  WAVE  DISPLACEMENT  SPECTRA 

We  consider  a  "large"  ocean  of  rectangular  area  A,  (periodic 
boundary  conditions)  with  a  plane  surface  at  z»0  (z>0  is  upl).  The 
vertical  displacement  from  equilibritun  is  C(3t, t).  This  is  expressed 
as  a  Fourier  expansion, 

C(x.t)  -  1  I  [ajjei(k*x-Wt)_c.c.]  .  (1.1) 


Here  x  -  (x,y)  and 

W(k)  -  [k(g+Yk2)]ii  (1.2) 

is  the  linear  wave  dispersion  relation.  We  use  mks  units,  so 

g»9.8m/s2  and  Y(*‘C/p)  is  7.5  x  10“5m^/s2  (t«7.5  x  10“2n/m  represents 

i 

a  nominal  value  for  the  surface  tension  of  uncontaroinated  water). 

For  linear  waves  the  a}^(t)  are  constant  in  time. 

I  The  point  of  writing  (1.1)  as  done  is  that  k  is  in  the  direc¬ 

tion  of  wave  propogation. 

I 

The  power  spectrum  or  displacement,  with  direction  of  propoga- 
I  tion  accounted  for,  is 

'  vKk)  -  [i<  I  a|.2  I  >]  (1.3) 

(W^  2 
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<C^>  ■  /d^k'Kk)  -  ^  i  <  I  aj  2 


>  . 


(1.4) 


The  symbol  represents  an  average  over  an  ensemble  of  oceans. 

The  spectrum  ({>  is  often  expressed  in  the  form 

vKk)  -  s(k)G(k.0),  (1.5) 

where  k  -  (k,0)  with  0-0  corresponding  to  the  wind  direction.  By 
convention, 

IT 

/Gd0  -  1  .  (1.6) 

-TT 

The  spectrum  is  also  expressed  in  terms  of  frequency  using  (1.2)  and 
the  relation 


S(k)  k  dk  -  Sf(«)  dW  .  (1.7) 

For  linear  waves  the  energy/unit  area  is 

E(k)  -  p“i  vKk)  ,  (1.8) 

k 

and  the  action  density  is 

F(k)  -  E(k)/«  .  (1.9) 
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To  describe  the  evolution  of  the  action  spectrum  F  in  the  pre¬ 
sence  of  a  large  scale  surface  current  Ug(x,t)  the  radiative 
transport  equation  is  frequently  used: 

[—  +  X  •  ^  -  Sjjj  +  Sy,  +  ,  (1.10) 

at 

Here  the  ray  trajectories  are  calculated  from  the  equations 

X  - 

k  -  -VjjH 

H  -  W(k)  +  k  •  Ug  .  (1.11) 

On  the  right-hand  side  of  (1.10),  represents  the  decay  rate 
due  to  viscosity,  Sy,  the  wind  growth  rate,  and  Sj.n  the  effects  of 
nonlinear  wave-wave  interaction. 
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2.0  GROWTH  AND  DECAY  OF  SURFACE  WAVES 

In  this  section  we  discuss  the  right-hand  side  of  (1.10).  The 
decay  rate  due  to  viscosity  is  [see,  for  example,  Phillips,  1977] 

Sv  -  -4Vk2F  .  (2.1) 

We  shall  use  here  the  nominal  value  of  1.1  x  10~^m^/s  for  the  kine¬ 
matic  viscosity  V.  Surface  contaminants  may  require  modification  of 
this  value. 

The  growth  due  to  the  wind  stress  is  described  by  the  term  S^,. 

The  wind  stress  on  the  surface  is 

Xw  -  Pa  (2.2) 

where  Pa  is  the  density  of  air  and  u*  is  the  friction  velocity.  We 
shall  use  here  the  values  given  by  Garratt  (1977): 

u*  -  U(10)[l  +  0.089  U(10)]^  »  (2.3) 

where  U(10)  is  the  wind  velocity  at  lOra  above  the  surface.  At  a 
height  z  we  shall  use  the  logorithroic  scaling  relation 

U(z)  -  ^ln(z/z,)  .  (2.4) 

K 
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This  is  appropriate  in  the  atmospheric  surface  boundary  layer  for 
conditions  of  neutral  stability  (a  condition  that  tends  to  be  valid 
over  oceans;  see  Panofsky  and  Dutton,  1984,  for  a  detailed 
discussion  and  references).  An  illustration  of  wind  flow  over  a 
Minnesota  wheat  field  (Kalmal,  et  al.  1976)  is  shown  in  Figure  2-1, 
where  the  planetary  boundary  layer  height  Zj  ■  1250m. 

The  quantity  K  a  0.4  in  (2.4)  is  the  von  Karman  constant  and 
Z.  is  the  "surface  roughness" 

Z.  -  0.014^*^  (2.5) 

g 

(Garratt,  1977).  This  value  is  appropriate  for  the  oceans.  Over 
land,  terrain  topology  can  lead  to  very  different  values  (see  for 
example,  Panofsky  and  Dutton,  1984). 

Several  recent  analyses  of  the  wave  growth  data  have  used  the 
form  implied  by  the  Miles  Theory, 

Sv  -  e(k)F.  .  (2.6) 

Empirical  models  for  6  are  deduced.  To  understand  the  basis  for 
doing  this,  we  briefly  review  the  Miles  Theory.  Write 
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c  -  eik-x  ^ 


^  -  eik-x  (>, 


(2.7) 


where  ^  is  the  velocity  potential.  The  linearized  Bernoulli 
equation  and  the  kinematic  boundary  condition  read 


+  I  +  (g  +  Vk2)C  -  0  , 


(2.8) 


A  linear  relation  is  postulated  to  relate  the  pressure  variation  to 
the  displacement, 

p,-  p(a  +  iPg)  (w2k)^  .  (2.9) 


On  replacing  ®  by  (-ifl),  we  obtain  the  dispersion  relation 
ot 

n  a  Wtl  +  i  (a  +  iPe)]  . 

2  * 


(2.10) 


The  rate  3  in  (2.6)  is  see.*,  to  be 
0  -  MgW  . 


(2.11) 
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Experimental  techniques  to  measure  wave  growth  vary,  and 
include  wave  tanks  and  field  measurements.  The  growth  may  be 
observed  directly,  using  wave  staffs,  laser  slope  meters,  electro¬ 
magnetic  waves,  etc.  (Donelan  et  al.  1985,  and,  Larson  and  Wright, 
1975).  An  alternative  method  is  to  measure  the  pressure  variation  p 
and  displacement  C*  Fourier  transform  these,  and  deduce  the  parame¬ 
ters  (a,  fig)  in  (2.9).  This  has  been  done,  for  example,  by  Snyder, 
et  al.  (1981),  Hsiao  and  Shemdin  (1983),  and  Hasselmann  et  al  (1986). 

An  extensive  review  of  wave  growth  data  published  prior  to  1980 
was  made  by  Plant  (1981).  This  data  included  wind  speeds  to  15m/s 
and  frequencies  in  the  range 

_ g _  <  W  <  40ir  .  (2.12) 

U(10) 

He  suggests  the  model 

0D  -  0.04(^)2«  cos  ♦  .  (2.13) 

Here  V  -  U/k  and  0  is  the  angle  between  wind  and  wave  direction. 
Evidence  for  the  factor  "cos  is  very  weak.  For  higher  wind 
speeds,  Amorocho  and  de  Vries  (1980)  describe  some  growth  rate  data. 
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Hsiao  and  Shemdin  (1983)  deduce  the  model 


Pg  -  1.4  X  10-^  (P-1)2  ‘  (2.14) 

]l  -  Q.aS[U(10)]  cos  4  .  (2.14) 

V 

This  is  based  on  measurements  in  the  range 

1  <  P  <  7,  5  <  U(10)  <  14m/s  .  (2.15) 

Donelan  and  Pierson  (1987)  proposed  a  model  valid  for  the  capillary 
range,  baaed  on  the  data  of  Larson  and  Wright  (1975).  This  tank 
data  included  wavelengths  A  in  the  range 

0.7  <  A  <  7.0cm  ,  (2.16) 


and 


0.17  <  ua  <  1 . 2ra/ s  . 


(2.16) 


Donelan  and  Pierson  purpose 


0  - 

D 


2.3  X 


10-4  [U(V2)  _  l]2 
V 


(2.17) 


The  use  of  tank  data  to  deduce  growth  rates  on  the  generally  much 
rougher  ocean  surface  has  an  unknown  validity. 
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A  comparison  of  the  above  growth  rates  is  illustrated  in 
Figures  (2),  (3),  and  (4).  The  principal  discrepency  is  for  A  in 
the  centimeter  range  for  which  (2.13)  and  (2.14)  are  not  expected  to 
be  valid.  For  long  wavelengths  the  agreement  is  fair.  For  waves  in 
the  10 's  of  meter  range  generation  by  shorter  waves  is  considered  by 
some  to  be  important. 

The  term  in  (1.10)  describes  the  effects  of  non-linear 
wave-wave  Interactions.  The  most  elaborate  model  for  this  has  been 
given  by  Hasselmann,  who  used  an  assumption  of  weak  interactions  and 
curomulant  closure  to  obtain  a  Boltzmann-like  integral.  Some  recent 
calculation  using  the  Hasselmann  theory  have  recently  been  published 
by  van  Gastel  (1987  a,b). 

The  complexity  and  uncertain  accuracy  of  the  Hasselmann  theory 
has  led  to  some  simplified  models.  A  very  simple  phenomenological 
model  was  suggested  by  Hughes  (1976),  which  was  generalized  by 
Phillips  (1985).  When  F  is  sufficiently  close  to  an  equilibrium 
form,  Fgq,  these  models  may  be  approximated  as 

St  -  -BiCF-Feq)  .  (2.18) 

Estimates  for  3t  were  made  by  Watson  (1986)  using  non-linear  wave 
theory.  These  are  reproduced  in  Figure  2-5. 
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Wave  Qrowth  Rate  (s-l) 


Figure  2-3.  Comparison  of  growth  rate  models  for  /3. 


Rates  (s-1) 


2- 


3.0  EQUILIBRIUM  SPECTRAL  MODELS 

Some  time  ago  Phillips  suggested  that  under  conditions  of  ade¬ 
quate  fetch  and  wind  duration  wave  spectra  tend  toward  an 
equilibrivim  state.  He  argued  that  there  should  be  an  equilibrixim 
range  between  waves  moving  at  wind  speed  [U(  10)i>V(k)  ]  and  the  region 
of  viscious  decay.  In  this  domain  he  proposed  that 

3  _  constant 

the  constant  being  dimensionless  (no  relevant  parameters!)  and  the 
power  of  k  determined  by  dimensional  arguments.  Kitaigorodskii  (in 
Phillips  and  Hasselmann,  1986)  has  recently  reviewed  the  philosophy 
of  equilibrium  spectrum  models. 

Pierson  and  Moskowitz  (1964)  proposed  a  more  elaborate  spectrum 
based  on  Phillips'  ideas: 

S  -  Sp(k)  -  (^JLiOli)  exp  [-.74(k*/k)2]  ,  (3.1) 

k^ 

where 

k*  -  _B -  .  (3.2) 

u2(10) 
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The  JONSWAP  experiment  led  to  the  replacement  of  the 


Pierson-Moskowitz  exponential  in  (3.1)  by  the  "peaked"  exponential 
e~^»  where 

r  *  0.74  (^*)^  -  0.5exp[-  “  0 . 9jAc*) ^ ]  .  (3.3) 

k  0.4kA 

Increasing  evidence  for  change  led  Phillips  (1986),  Donelan  et 
al.  (1985),  and  others  to  give  up  the  k~^  spectral  form  by  a  factor  of 
J^.  They  proposed  that 


S  -  SD(k) 


_ a _ 

a-7/2k*-l/2  e-r 


(3.4) 


where  a  reasonable  choice  for  the  dimensional  constant  *  is 


a  -  2  X  10-3  , 


(3.5) 


Donalen  et  al.  suggest  some  generalization  of  (3.4)  for  limited  fetch 
conditions . 

Observation  by  M.  Banner  (private  communication,  see  also 
Donalen  and  Pierson,  1987)  suggest  that  for  k^  <  k  <  k2,  a  reaso¬ 
nable  model  for  S  is 

S  -  SB(k)  -  ,  (3.6) 

k^ 
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where  a 'is  chosen  to  give 

SD(ki)  -  SB(ki) 


(3.7) 


at 


ki  ■  4m  ^  . 

A  number  of  observations  suggest  that  for 


(3.7) 


k  >  k2  “  200m” ^ 


(3.8) 


further  models  are  needed  for  S.  Bjerkaas  and  Riedel  (1979)  have 
reviewed  the  data  (particularly  that  of  Mitsuyasu,  1977)  and  have 
developed  an  elaborate  model  for  the  range  (3.8).  A  simplified  ver¬ 
sion  of  this  model  is 


S  -  SR(k) 


4.4  X  10-5  (2^“^ 
k[k(l+^-^)  (l+exp(k-125u*)  ] 


3k2 


(3.9) 


k  >  k2  • 


(3.9) 


Here 


p  »  3  -  0.4341n(u*) 


(3.10) 
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with  ua  in  m/s  and 


The  resulting  spectrxim  using  Sj),  Sg,  and  Sg  in  the  ranges 
described  is  illustrated  in  Figures  (6),  (7),  and  (8).  Donelan  and 
Pierson  (1987)  model  the  regime  k  >  k2  without  using  (3.9),  but  a 
version  of  (3.6). 

Models  for  the  approach  to  equilibrium  have  been  suggested  by 
Hasselmann  et  al  (1976).  An  example  is  given  for  which 

U(10)  -  0,  t  <  0 

»  U,,  a  constant  for  t  >  0  .  (3.12) 


A  parameter  is  defined 


3 

-  120  (&t)'7  ,  >  1  , 

U 


-  1,  if  above  is  less  than  unity. 


(3.13) 


Then  T  in  (3.4)  is  modified  by  replacing  k*  by 


kpj  “  ' 


(3.14) 
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>1.0001 


Wind  speed  =  4m/s 


Figure  3*1.  Composite  spectral  model  for  S(k). 
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Ms. 0001 


20 


Wind  speed  =  7  m/s 


Figure  3-2.  Composite  spectral  model  for  S(k). 

I 
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1000*k**4*S(k) 


Wind  speed  =•  lOm/s 


Figure  3-3.  Composite  spectral  model  for  S(k). 
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A  more  recent  review  has  led  Donelan  et  al.  (1985)  to  suggest 
replacing  (3.16)  by 


G  -  2  sech2(O0)  ,  (3.18) 

2 

2 

0  -  2.9(^*)3  ,  1  <  ^  <  4  , 

TT  k* 
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1.2,  otherwise. 


(3.19) 


(Equation  (3.19)  represents  my  simplification  of  a  more  elaborate 
representation. ) 

The  diversity  of  spectral  models  and  the  recent  dates  on  many 
of  the  references  will  convince  you  that  further  models  can  be 
expected.  Comparison  of  the  models  described  above  suggests  that 
changes  have  tended  to  be  more  evolutionary  than  revolutionary  in 
this  field  and  that  the  existing  models  can  be  useful  even  if  of 
limited  precision. 
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